Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico (preprint)

Up to November 2021, over 200 different SARS-CoV-2 lineages circulated in Mexico. To investigate lineage replacement dynamics, we applied a phylodynamic approach to explore the evolutionary trajectories of five dominant lineages that circulated during the first year of the local epidemic. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in the country. Lineages B.1.1.222 and B.1.1.519 showed comparable dynamics, represented by clades likely originating in Mexico and persisting for over a year. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. We further explored viral movements across the country, applied within the largest clades identified (belonging to lineage B.1.617.2). Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

Comparing the evolutionary dynamics of predominant SARS-CoV-2 virus lineages co-circulating in Mexico (preprint)

Up to November 2021, over 200 different SARS-CoV-2 lineages circulated in Mexico. To investigate lineage replacement dynamics, we applied a phylodynamic approach to explore the evolutionary trajectories of five dominant lineages that circulated during the first year of the local epidemic. For most lineages, peaks in sampling frequencies coincided with different epidemiological waves of infection in the country. Lineages B.1.1.222 and B.1.1.519 showed comparable dynamics, represented by clades likely originating in Mexico and persisting for over a year. Lineages B.1.1.7, P.1 and B.1.617.2 also displayed similar dynamics, characterized by multiple introduction events leading to a few successful extended local transmission chains that persisted for several months. We further explored viral movements across the country, applied within the largest clades identified (belonging to lineage B.1.617.2). Many clades were located within the south region of the country, suggesting that this area played a key role in the spread of SARS-CoV-2 in Mexico.

The origins and molecular evolution of SARS-CoV-2 lineage B.1.1.7 in the UK (preprint)

The first SARS-CoV-2 variant of concern (VOC) to be designated was lineage B.1.1.7, later labelled by the World Health Organisation (WHO) as Alpha. Originating in early Autumn but discovered in December 2020, it spread rapidly and caused large waves of infections worldwide. The Alpha variant is notable for being defined by a long ancestral phylogenetic branch with an increased evolutionary rate, along which only two sequences have been sampled. Alpha genomes comprise a well-supported monophyletic clade within which the evolutionary rate is more typical of SARS-CoV-2. The Alpha epidemic continued to grow despite the continued restrictions on social mixing across the UK, and the imposition of new restrictions, in particular the English national lockdown in November 2020. While from a case-number perspective these interventions succeeded in reducing the absolute number of cases of SARS-CoV-2 in the UK, the impact of these non-pharmaceutical interventions was predominantly to drive the decline of those SARS-CoV-2 lineages that preceded Alpha. We investigate the only two sampled sequences that fall on the branch ancestral to Alpha. We find that one is likely to be a true intermediate sequence, providing information about the order of mutational events that led to Alpha. We explore alternate hypotheses that can explain how Alpha acquired a large number of mutations yet remained largely unobserved in a region of high genomic surveillance: an under-sampled geographical location, a non-human animal population, or a chronically-infected individual. We conclude that the last hypothesis provides the best explanation of the observed behaviour and dynamics of the variant, although we find that the individual need not be immunocompromised, as persistently-infected immunocompetent hosts also display a higher within-host rate of evolution. Finally, we compare the ancestral branches and mutation profiles of other VOCs to each other, and identify that Delta appears to be an outlier both in terms of the genomic locations of its defining mutations, and its lack of rapid evolutionary rate on the ancestral branch. As new variants, such as Omicron, continue to evolve (potentially through similar mechanisms) it remains important to investigate the origins of other variants to identify ways to potentially disrupt their evolution and emergence.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant (preprint)

The Delta variant of concern of SARS-CoV-2 has spread globally causing large outbreaks and resurgences of COVID-19 cases. The emergence of Delta in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 Delta genomes from England in combination with 93,649 global genomes to reconstruct the emergence of Delta, and quantify its introduction to and regional dissemination across England, in the context of changing travel and social restrictions. Through analysis of human movement, contact tracing, and virus genomic data, we find that the focus of geographic expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced >1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers from India reduced onward transmission from importations; however the transmission chains that later dominated the Delta wave in England had been already seeded before restrictions were introduced. In England, increasing inter-regional travel drove Delta's nationwide dissemination, with some cities receiving >2,000 observable lineage introductions from other regions. Subsequently, increased levels of local population mixing, not the number of importations, was associated with faster relative growth of Delta. Among US states, we find that regions that previously experienced large waves also had faster Delta growth rates, and a model including interactions between immunity and human behaviour could accurately predict the rise of Delta there. Delta’s invasion dynamics depended on fine scale spatial heterogeneity in immunity and contact patterns and our findings will inform optimal spatial interventions to reduce transmission of current and future VOCs such as Omicron.

Emergence and widespread circulation of a recombinant SARS-CoV-2 lineage in North America (preprint)

Genetic recombination is an important driving force of coronavirus evolution. While some degree of virus recombination has been reported during the COVID-19 pandemic, previously detected recombinant lineages of SARS-CoV-2 have shown limited circulation and been observed only in restricted areas. Prompted by reports of unusual genetic similarities among several Pango lineages detected mainly in North and Central America, we present a detailed phylogenetic analysis of four SARS-CoV-2 lineages (B.1.627, B.1.628, B.1.631 and B.1.634) in order to investigate the possibility of virus recombination among them. Two of these lineages, B.1.628 and B.1.631, are split into two distinct clusters (here named major and minor). Our phylogenetic and recombination analyses of these lineages find well-supported phylogenetic differences between the Orf1ab region and the rest of the genome (S protein and remaining reading frames). The lineages also contain several deletions in the NSP6, Orf3a and S proteins that can augment reconstruction of reliable evolutionary histories. By reconciling the deletions and phylogenetic data, we conclude that the B.1.628 major cluster originated from a recombination event between a B.1.631 major virus and a lineage B.1.634 virus. This scenario inferred from genetic data is supported by the spatial and temporal distribution of the three lineages, which all co-circulated in the USA and Mexico during 2021, suggesting this region is where the recombination event took place. We therefore support the designation of the B.1.628 major cluster as recombinant lineage XB in the Pango nomenclature. The widespread circulation of lineage XB across multiple countries over a longer timespan than the previously designated recombinant XA lineage raises important questions regarding the role and potential effects of recombination on the evolution of SARS-CoV-2 during the ongoing COVID-19 pandemic.

Genomic epidemiology of early SARS-CoV-2 transmission dynamics in Gujarat, India (preprint)

Genomic surveillance of SARS-CoV-2 has played a decisive role in understanding the transmission and evolution of the virus during its emergence and continued circulation. However, limited genomic sampling in many high-incidence countries has impeded detailed studies of SARS-CoV-2 genomic epidemiology. Consequently, critical questions remain about the generation and global distribution of virus genetic diversity. To address this gap, we investigated SARS-CoV-2 transmission dynamics in Gujarat, India, during its first epidemic wave and shed light on virus spread in one of the pandemics hardest-hit regions. By integrating regional case data and 434 whole virus genome sequences sampled across 20 districts from March to July 2020, we reconstructed the epidemic dynamics and spatial spread of SARS-CoV-2 in Gujarat, India. Our findings revealed that global and regional connectivity, along with population density, were significant drivers of the Gujarat SARS-CoV-2 outbreak. The three most populous districts in Gujarat accounted [~]84% of total cases during the first wave. Moreover, we detected over 100 virus lineage introductions, which were primarily associated with international travel. Within Gujarat, virus dissemination occurred predominantly from densely populated regions to geographically proximate locations with low-population density. Our findings suggest SARS-CoV-2 transmission follows a gravity model in India, with urban centres contributing disproportionately to onward virus spread.

Global disparities in SARS-CoV-2 genomic surveillance (preprint)

Genomic sequencing provides critical information to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments and vaccines, and guide public health responses. To investigate the spatiotemporal heterogeneity in the global SARS-CoV-2 genomic surveillance, we estimated the impact of sequencing intensity and turnaround times (TAT) on variant detection in 167 countries. Most countries submit genomes >21 days after sample collection, and 77% of low and middle income countries sequenced <0.5% of their cases. We found that sequencing at least 0.5% of the cases, with a TAT <21 days, could be a benchmark for SARS-CoV-2 genomic surveillance efforts. Socioeconomic inequalities substantially impact our ability to quickly detect SARS-CoV-2 variants, and undermine the global pandemic preparedness. One-Sentence SummarySocioeconomic inequalities impacted the SARS-CoV-2 genomic surveillance, and undermined the global pandemic preparedness.

Pango lineage designation and assignment using SARS-CoV-2 spike gene nucleotide sequences (preprint)

More than 2 million SARS-CoV-2 genome sequences have been generated and shared since the start of the COVID-19 pandemic and constitute a vital information source that informs outbreak control, disease surveillance, and public health policy. The Pango dynamic nomenclature is a popular system for classifying and naming genetically-distinct lineages of SARS-CoV-2, including variants of concern, and is based on the analysis of complete or near-complete virus genomes. However, for several reasons, nucleotide sequences may be generated that cover only the spike gene of SARS-CoV-2. It is therefore important to understand how much information about Pango lineage status is contained in spike-only nucleotide sequences. Here we explore how Pango lineages might be reliably designated and assigned to spike-only nucleotide sequences. We survey the genetic diversity of such sequences, and investigate the information they contain about Pango lineage status. Although many lineages, including the main variants of concern, can be identified clearly using spike-only sequences, some spike-only sequences are shared among tens or hundreds of Pango lineages. To facilitate the classification of SARS-CoV-2 lineages using subgenomic sequences we introduce the notion of designating such sequences to a "lineage set", which represents the range of Pango lineages that are consistent with the observed mutations in a given spike sequence. These data provide a foundation for the development of software tools that can assign newly-generated spike nucleotide sequences to Pango lineage sets.

Purifying selection determines the short-term time dependency of evolutionary rates in SARS-CoV-2 and pH1N1 influenza (preprint)

High throughput sequencing enables rapid genome sequencing during infectious disease outbreaks and provides an opportunity to quantify the evolutionary dynamics of pathogens in near real-time. One difficulty of undertaking evolutionary analyses over short timescales is the dependency of the inferred evolutionary parameters on the timespan of observation. Here, we characterise the molecular evolutionary dynamics of SARS-CoV-2 and 2009 pandemic H1N1 (pH1N1) influenza during the first 12 months of their respective pandemics. We use Bayesian phylogenetic methods to estimate the dates of emergence, evolutionary rates, and growth rates of SARS-CoV-2 and pH1N1 over time and investigate how varying sampling window and dataset sizes affects the accuracy of parameter estimation. We further use a generalised McDonald-Kreitman test to estimate the number of segregating non-neutral sites over time. We find that the inferred evolutionary parameters for both pandemics are time-dependent, and that the inferred rates of SARS-CoV-2 and pH1N1 decline by ~50% and ~100%, respectively, over the course of one year. After at least 4 months since the start of sequence sampling, inferred growth rates and emergence dates remain relatively stable and can be inferred reliably using a logistic growth coalescent model. We show that the time-dependency of the mean substitution rate is due to elevated substitution rates at terminal branches which are 2-4 times higher than those of internal branches for both viruses. The elevated rate at terminal branches is strongly correlated with an increasing number of segregating non-neutral sites, demonstrating the role of purifying selection in generating the time-dependency of evolutionary parameters during pandemics.

Viral infection and transmission in a large, well-traced outbreak caused by the SARS-CoV-2 Delta variant (preprint)

We report the first local transmission of the SARS-CoV-2 Delta variant in mainland China. All 167 infections could be traced back to the first index case. Daily sequential PCR testing of the quarantined subjects indicated that the viral loads of Delta infections, when they first become PCR+, were on average ~1000 times greater compared to A/B lineage infections during initial epidemic wave in China in early 2020, suggesting potentially faster viral replication and greater infectiousness of Delta during early infection. We performed high-quality sequencing on samples from 126 individuals. Reliable epidemiological data meant that, for 111 transmission events, the donor and recipient cases were known. The estimated transmission bottleneck size was 1-3 virions with most minor intra-host single nucleotide variants (iSNVs) failing to transmit to the recipients. However, transmission heterogeneity of SARS-CoV-2 was also observed. The transmission of minor iSNVs resulted in at least 4 of the 30 substitutions identified in the outbreak, highlighting the contribution of intra-host variants to population level viral diversity during rapid spread. Disease control activities, such as the frequency of population testing, quarantine during pre-symptomatic infection, and level of virus genomic surveillance should be adjusted in order to account for the increasing prevalence of the Delta variant worldwide.

Viral infection and transmission in a large, well-traced outbreak caused by the SARS-CoV-2 Delta variant (preprint)

Summary We report the first local transmission of the SARS-CoV-2 Delta variant in mainland China. All 167 infections could be traced back to the first index case. Daily sequential PCR testing of the quarantined subjects indicated that the viral loads of Delta infections, when they first become PCR+, were on average ∼1000 times greater compared to A/B lineage infections during initial epidemic wave in China in early 2020, suggesting potentially faster viral replication and greater infectiousness of Delta during early infection. We performed high-quality sequencing on samples from 126 individuals. Reliable epidemiological data meant that, for 111 transmission events, the donor and recipient cases were known. The estimated transmission bottleneck size was 1-3 virions with most minor intra-host single nucleotide variants (iSNVs) failing to transmit to the recipients. However, transmission heterogeneity of SARS-CoV-2 was also observed. The transmission of minor iSNVs resulted in at least 4 of the 30 substitutions identified in the outbreak, highlighting the contribution of intra-host variants to population level viral diversity during rapid spread. Disease control activities, such as the frequency of population testing, quarantine during pre-symptomatic infection, and level of virus genomic surveillance should be adjusted in order to account for the increasing prevalence of the Delta variant worldwide.

Genomics and epidemiology of a novel SARS-CoV-2 lineage in Manaus, Brazil (preprint)

Cases of SARS-CoV-2 infection in Manaus, Brazil, resurged in late 2020, despite high levels of previous infection there. Through genome sequencing of viruses sampled in Manaus between November 2020 and January 2021, we identified the emergence and circulation of a novel SARS-CoV-2 variant of concern, lineage P.1, that acquired 17 mutations, including a trio in the spike protein (K417T, E484K and N501Y) associated with increased binding to the human ACE2 receptor. Molecular clock analysis shows that P.1 emergence occurred around early November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.4–2.2 times more transmissible and 25-61% more likely to evade protective immunity elicited by previous infection with non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness. One-Sentence Summary We report the evolution and emergence of a SARS-CoV-2 lineage of concern associated with rapid transmission in Manaus.

Levels of SARS-CoV-2 Lineage P.1 Neutralization by Antibodies Elicited after Natural Infection and Vaccination (preprint)

Background: A new SARS-CoV-2 lineage, named P.1 (20J/501Y.V3), has recently been detected in Brazil. Mutations accrued by the P.1 lineage include amino acid changes in the receptor-binding domain of the spike protein that also are reported in variants of concern in the United Kingdom (B.1.1.7) and South Africa (B.1.325).Methods: We isolated two P.1-containing specimens from nasopharyngeal and bronchoalveolar lavage samples of patients of Manaus, Brazil. We measured neutralization of the P.1 virus after incubation with the plasma of 19 COVID-19 convalescent blood donors and recipients of the chemically-inactivated CoronaVac vaccine and compared these results to neutralization of a SARS-CoV-2 B-lineage previously circulating in Brazil.Findings: The immune plasma of COVID-19 convalescent blood donors had 6-fold less neutralizing capacity against the P.1 than against the B-lineage. Moreover, five months after booster immunization with CoronaVac, plasma from vaccinated individuals failed to efficiently neutralize P.1 lineage isolates.Interpretation: These data indicate that the P.1 lineage may escape from neutralizing antibodies generated in response to polyclonal stimulation against previously circulating variants of SARS-CoV-2.Funding: São Paulo Research Foundation, MCTI/FINEP, Medical Research Council, National Council for Scientific and Technological Development, National Institutes of Health.Conflict of Interest: M.S.D. is a consultant for Inbios, Vir Biotechnology, NGMBiopharmaceuticals, and Carnival Corporation, and on the Scientific Advisory Boards of Moderna and Immunome. The Diamond laboratory has received funding support in sponsored research agreements from Moderna, Vir Biotechnology, and Emergent BioSolutions.Ethical Approval: All procedures followed the ethical standards of the responsible committee on humanexperimentation and approved by the ethics committees from the University of Campinas, Brazil (Approval number CONEP 4.021.484 for plasma collection of blood donors, CAEE32078620.4.0000.5404 and 30227920.9.0000.5404 for the sampling of vaccinated and viral genome sequencing, respectively). All patient data were anonymized before study inclusion.Informed consent was obtained from all subjects for being included in the study.

Epidemiological and evolutionary considerations of SARS-CoV-2 vaccine dosing regimes (preprint)

As the threat of Covid-19 continues and in the face of vaccine dose shortages and logistical challenges, various deployment strategies are being proposed to increase population immunity levels. How timing of delivery of the second dose affects infection burden but also prospects for the evolution of viral immune escape are critical questions. Both hinge on the strength and duration (i.e. robustness) of the immune response elicited by a single dose, compared to natural and two-dose immunity. Building on an existing immuno-epidemiological model, we find that in the short-term, focusing on one dose generally decreases infections, but longer-term outcomes depend on this relative immune robustness. We then explore three scenarios of selection, evaluating how different second dose delays might drive immune escape via a build-up of partially immune individuals. Under certain scenarios, we find that a one-dose policy may increase the potential for antigenic evolution. We highlight the critical need to test viral loads and quantify immune responses after one vaccine dose, and to ramp up vaccination efforts throughout the world.

Establishment & lineage dynamics of the SARS-CoV-2 epidemic in the UK (preprint)

The UK's COVID-19 epidemic during early 2020 was one of world's largest and unusually well represented by virus genomic sampling. Here we reveal the fine-scale genetic lineage structure of this epidemic through analysis of 50,887 SARS-CoV-2 genomes, including 26,181 from the UK sampled throughout the country's first wave of infection. Using large-scale phylogenetic analyses, combined with epidemiological and travel data, we quantify the size, spatio-temporal origins and persistence of genetically-distinct UK transmission lineages. Rapid fluctuations in virus importation rates resulted in >1000 lineages; those introduced prior to national lockdown were larger and more dispersed. Lineage importation and regional lineage diversity declined after lockdown, whilst lineage elimination was size-dependent. We discuss the implications of our genetic perspective on transmission dynamics for COVID-19 epidemiology and control.

Clinical, Immunological and Virological Characterization of COVID-19 Patients that Test Re-Positive for SARS-CoV-2 by RT-PCR (preprint)

Background: COVID-19 pandemic is underway. Some COVID-19 cases re-tested positive for SARS-CoV-2 RNA after discharge raising the public concern on their infectivity. Characterization of re-positive cases are urgently needed for designing intervention strategies. Methods: Clinical data were obtained through Guangdong COVID-19 surveillance network. Neutralization antibody titre was determined using a microneutralization assay. Potential infectivity of clinical samples was evaluated after the cell inoculation. SARS-CoV-2 RNA was detected using three different RT-PCR kits and multiplex PCR with nanopore sequencing. Findings: Among 619 discharged COVID-19 cases, 87 were re-tested as SARS-CoV-2 positive in circumstance of social isolation. All re-positive cases had mild or moderate symptoms in initial diagnosis and a younger age distribution (mean, 30·4). Re-positive cases (n=59) exhibited similar neutralization antibodies (NAbs) titre distributions to other COVID-19 cases (n=150) parallel-tested in this study. No infective viral strain could be obtained by culture and none full-length viral genomes could be sequenced for all re-positive cases. Interpretation: Re-positive SARS-CoV-2 was not caused by the secondary infection and was identified in around 14% of discharged cases. A robust Nabs response and a potential virus genome degradation were detected from nearly all re-positive cases suggesting a lower transmission risk, especially through a respiratory route. Funding: This work was supported by grants from Guangdong Provincial Novel Coronavirus Scientific and Technological Project (2020111107001), Science and Technology Planning Project of Guangdong(2018B020207006), National Science and Technology Project(2020YFC0846800).Declaration of Interests: All authors: No reported conflicts of interest.Ethics Approval Statement: This study was reviewed and approved by the Medical Ethical Committee of Guangdong Provincial Center for Disease Control and Prevention. Data collection and analysis of cases were determined by the Health Commission of Guangdong province to be part of a continuing public health outbreak investigation during the emergency response and were thus considered exempt from institutional review board approval.

Evaluating the effects of SARS-CoV-2 Spike mutation D614G on transmissibility and pathogenicity (preprint)

Global dispersal and increasing frequency of the SARS-CoV-2 Spike protein variant D614G are suggestive of a selective advantage but may also be due to a random founder effect. We investigate the hypothesis for positive selection of Spike D614G in the United Kingdom using more than 25,000 whole genome SARS-CoV-2 sequences. Despite the availability of a large data set, well represented by both Spike 614 variants, not all approaches showed a conclusive signal of positive selection. Population genetic analysis indicates that 614G increases in frequency relative to 614D in a manner consistent with a selective advantage. We do not find any indication that patients infected with the Spike 614G variant have higher COVID-19 mortality or clinical severity, but 614G is associated with higher viral load and younger age of patients. Significant differences in growth and size of 614G phylogenetic clusters indicate a need for continued study of this variant.

Clinical, immunological and virological characterization of COVID-19 patients that test re-positive for SARS-CoV-2 by RT-PCR (preprint)

BackgroundCOVID-19 pandemic is underway. Some COVID-19 cases re-tested positive for SARS-CoV-2 RNA after discharge raising the public concern on their infectivity. Characterization of re-positive cases are urgently needed for designing intervention strategies. MethodsClinical data were obtained through Guangdong COVID-19 surveillance network. Neutralization antibody titre was determined using a microneutralization assay. Potential infectivity of clinical samples was evaluated after the cell inoculation. SARS-CoV-2 RNA was detected using three different RT-PCR kits and multiplex PCR with nanopore sequencing. ResultsAmong 619 discharged COVID-19 cases, 87 were re-tested as SARS-CoV-2 positive in circumstance of social isolation. All re-positive cases had mild or moderate symptoms in initial diagnosis and a younger age distribution (mean, 30.4). Re-positive cases (n=59) exhibited similar neutralization antibodies (NAbs) titre distributions to other COVID-19 cases (n=150) parallel-tested in this study. No infective viral strain could be obtained by culture and none full-length viral genomes could be sequenced for all re-positive cases. ConclusionsRe-positive SARS-CoV-2 was not caused by the secondary infection and was identified in around 14% of discharged cases. A robust Nabs response and a potential virus genome degradation were detected from nearly all re-positive cases suggesting a lower transmission risk, especially through a respiratory route.

Epidemiological and clinical characteristics of the early phase of the COVID-19 epidemic in Brazil (preprint)

BackgroundThe first case of COVID-19 was detected in Brazil on February 25, 2020. We report the epidemiological, demographic, and clinical findings for confirmed COVID-19 cases during the first month of the epidemic in Brazil. MethodsIndividual-level and aggregated COVID-19 data were analysed to investigate demographic profiles, socioeconomic drivers and age-sex structure of COVID-19 tested cases. Basic reproduction numbers (R0) were investigated for Sao Paulo and Rio de Janeiro. Multivariate logistic regression analyses were used to identify symptoms associated with confirmed cases and risk factors associated with hospitalization. Laboratory diagnosis for eight respiratory viruses were obtained for 2,429 cases. FindingsBy March 25, 1,468 confirmed cases were notified in Brazil, of whom 10% (147 of 1,468) were hospitalised. Of the cases acquired locally (77{middle dot}8%), two thirds (66{middle dot}9% of 5,746) were confirmed in private laboratories. Overall, positive association between higher per capita income and COVID-19 diagnosis was identified. The median age of detected cases was 39 years (IQR 30-53). The median R0 was 2{middle dot}9 for Sao Paulo and Rio de Janeiro. Cardiovascular disease/hypertension were associated with hospitalization. Co-circulation of six respiratory viruses, including influenza A and B and human rhinovirus was detected in low levels. InterpretationSocioeconomic disparity determines access to SARS-CoV-2 testing in Brazil. The lower median age of infection and hospitalization compared to other countries is expected due to a younger population structure. Enhanced surveillance of respiratory pathogens across socioeconomic statuses is essential to better understand and halt SARS-CoV-2 transmission. FundingSao Paulo Research Foundation, Medical Research Council, Wellcome Trust and Royal Society.

Crowding and the epidemic intensity of COVID-19 transmission (preprint)

The COVID-19 pandemic is straining public health systems worldwide and major non-pharmaceutical interventions have been implemented to slow its spread. During the initial phase of the outbreak the spread was primarily determined by human mobility. Yet empirical evidence on the effect of key geographic factors on local epidemic spread is lacking. We analyse highly-resolved spatial variables for cities in China together with case count data in order to investigate the role of climate, urbanization, and variation in interventions across China. Here we show that the epidemic intensity of COVID-19 is strongly shaped by crowding, such that epidemics in dense cities are more spread out through time, and denser cities have larger total incidence. Observed differences in epidemic intensity are well captured by a metapopulation model of COVID-19 that explicitly accounts for spatial hierarchies. Densely-populated cities worldwide may experience more prolonged epidemics. Whilst stringent interventions can shorten the time length of these local epidemics, although these may be difficult to implement in many affected settings.